load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

location_f="/Users/ivana/WORK_BSC/AS/files/"

fw    = addfile(location_f + "seaiceconc_ano_nov_ipsl.nc", "r")
fc    = addfile(location_f + "seaiceconc_ttest_nov_ipsl.nc", "r")

print(fw)
print(fc)

;*******************************************;
;Read in variables
;*******************************************
lat=fw->LAT
lon=fw->LON
;print(lat)


nlat  = dimsizes(lat)
nlon  = dimsizes(lon)

ano=fw->SICF_ANO(0,:,:)
tempsc=fc->SEUIL_SICF_NOV

ano=ano*100. ; convert to %
printVarSummary(ano)


;this is not necessary but makes the plot nicer
do i=0,nlat-1
do j=0,nlon-1
if (.not.ismissing(ano(i,j)) ) then
if (abs(ano(i,j)).lt.0.00001)
ano(i,j)=ano@_FillValue
end if
else
ano(i,j)=ano@_FillValue
end if

end do
end do


prob=tempsc

;-------------------------
; making the plots
;-------------------------
wks = gsn_open_wks("pdf","NHallanoice_95per_nov_medpolar")

gsn_define_colormap(wks,"BlueDarkOrange18")

  plot = new(1,graphic)                          ; create a plot array

  res          = True
  res@gsnDraw  = False                          ; don't draw
  res@gsnFrame = False                          ; don't advance frame, these 2 are going to be set with panel plot - gsn_panel!
  res@cnInfoLabelOn = False                     ; turn off cn info label
  res@cnLinesOn           = False
  res@cnFillOn            = True          ; turn on color
  res@lbLabelBarOn        = True           ; turn off individual cb's

  res@mpPerimOn = True
  res@cnLevelSelectionMode =  "ExplicitLevels"

 res@cnLineLabelsOn=False
 res@mpCenterLonF     = 300.

 res@mpFillOn = True
 res@mpLandFillColor = "grey"
 res@mpFillDrawOrder = "PostDraw"

 res@gsnPolar   = "NH"
 res@mpMinLatF            = 50.
 
  res1          = True
  res1@gsnDraw  = False                          ; don't draw
  res1@gsnFrame = False                          ; don't advance frame, these 2 are going to be set with panel plot - gsn_panel!
  res1@cnInfoLabelOn = False                     ; turn off cn info label
  res1@cnLinesOn           = False; True
  res1@cnFillOn            = True          ; turn on color
  res1@lbLabelBarOn        = False           ; turn off individual cb's

  res@cnMissingValFillColor   = "white"
  res1@cnMonoFillPattern = False ;True
  res1@cnFillPatterns      = (/-1,8,-1/) ;17, 8, 13
  res1@cnLevelSelectionMode =  "ExplicitLevels"
  res1@cnLevels =(/0.99, 1.05/) 
  res1@cnMonoFillColor     = True         ; Use same fill color

 res1@cnLineLabelsOn=False
 res1@cnMonoLineColor= True
 res1@cnLineColor      = (/"grey36"/) ; (/"white"/)

res@mpGeophysicalLineThicknessF = 2.
res@lbAutoManage = False
res@lbLabelFontHeightF  = 0.017
res@pmLabelBarHeightF  = 0.08


res1@gsnCenterString = ""
res1@gsnLeftString = ""
res1@gsnRightString = ""

res@gsnLeftString = ""
res@gsnCenterString = ""

res@cnFillColors = (/2,3,4,5,7,9,11,12,13,14,15,16/)
res@cnLevels =(/-60.,-40.,-20.,-10.,-5.,-2.5,0,2.5, 5./)

res@gsnRightString = ""
plot(0) = gsn_csm_contour_map(wks,ano(:,:),res)
plot1 = gsn_csm_contour(wks,gsn_add_cyclic_point(prob),res1)
overlay(plot(0),plot1)

;________________________
; creating the panel:
;________________________
 resP            = True

resP@gsnPanelLabelBar    = False;True                ; add common colorbar
resP@lbLabelFontHeightF  = 0.007               ; make labels smaller

resP@gsnPanelBottom   = 0.05                   ; add space at bottom
;resP@gsnPanelFigureStrings= (/"a)","b)","c)"/) ; add strings to panel
res@txFontHeightF     = .24


resP@gsnMaximize    = True                ; maximize plot
gsn_panel(wks,plot,(/1,1/),resP)

